Velocity oscillations in actin-based motility 
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We present a simple and generic theoretical description of actin-based motility, where polymer- 
ization of filaments maintains propulsion. The dynamics is driven by polymerization kinetics at the 
filaments' free ends, crosslinking of the act in network, attachment and detachment of filaments to 
the obstacle interfaces and entropic forces. We show that spontaneous oscillations in the velocity 
emerge in a broad range of parameter values, and compare our findings with experiments. 
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Force generation by semiflexible polymers is versatilely 
used for cell motility. The leading edge of lamellipodia 
of crawling cells pQ is pushed forward by a polymerizing 
actin network and bacteria move inside cells by riding on 
a comet tail of growing actin filaments [2 , 3 . In vivo sys- 
tems are complemented by in vitro assays using plastic 
beads and lipid vesicles [4[ [5j [6] . The defining feature of 
semiflexible polymers is the order of magnitude of their 
bending energy which is in the range of k^T. They un- 
dergo thermal shape fluctuations and the force exerted 
by the filaments against an obstacle arises from elastic 
and entropic contributions [7J [8] . 

Mathematical models have quantified the force gener- 
ated by actin filaments growing against obstacles [7] [9] . 
The resisting force depends on the obstacle which is 
pushed. In case of pathogens, it has a small com- 
ponent from viscous drag of the moving obstacle but 
consists mainly of the force exerted by actin filaments 
bound to the surface of the bacteria and pulling it back- 
wards p~0j[TT]. The tethered ratchet model [12 is a math- 
ematical formulation of these experimental findings in 
terms of the dynamics of the number of attached and de- 
tached polymers. The starting point of our approach will 
be the dynamics of the distributions of the free length of 
both polymer populations. 

Actin polymerization in living cells and extracts is con- 
trolled by a complex molecular network [3]. Nucleation 
of new filaments, capping of existing ones, exchange of 
ADP for ATP on actin monomers, buffering of monomers 
etc. all contribute to that control and have been modeled 
[12] [13] [14] . Our goal is not to model the full complex- 
ity of that biochemical network. Rather we focus on the 
core process of force generation and force balance ensuing 
from the interplay between bound pulling filaments and 
polymerizing pushing filaments, the transition between 
these two groups and the motion of the whole force gen- 
erating configuration. This is motivated by recent ob- 
servations of complex dynamics in simple reconstituted 
systems: the velocity of beads or pathogens propelled 
by actin polymerization may oscillate [16] [I7j [18] . Our 
goal is to describe the dynamics of such biochemically 
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FIG. 1: (color online) Schematic representation of an ensem- 
ble of actin filaments oriented at $ — with respect to the 
normal h of an obstacle interface, which may either be a cell 
membrane or a bacterium. While attached filaments are un- 
der tension and pull the interface back, detached filaments 
are compressed, elongate by polymerization with rate k on 
and push the interface forward. All filaments in the brush 
are firmly anchored in a cross- linked network, whose front 
advances with velocity v g reducing the free length / of the 
filaments. Attached filaments detach with stress dependent 
rate kd and detached filaments attach with constant rate k a - 
v is the interface velocity in the extracellular medium, and 
x is the distance between the front of the network and the 
interface. 



simpler systems and find a robust microscopic descrip- 
tion for oscillation mechanisms, which may then be con- 
trolled by higher order processes. Such a study is meant 
to complement investigations based on a continuum ap- 
proach pa eng. 

We consider a fixed number N of actin filaments [19 
firmly anchored into a rigid cross-linked network, which 
advances with velocity v g \ for an illustration see Fig. [I] 
Filaments of variable length / are either attached to the 
obstacle interface via a protein complex or detached from 
it, with time-dependent number distributions denoted by 
N a (l,t) and Nd(l,i), respectively. In the detached state, 
filaments polymerize at a velocity v v {l,x), which depends 
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on both the polymer length I and the distance x between 
rigid support and obstacle. Transitions between the two 
filament populations occur with a constant attachment 
rate k a and a stress-dependent detachment rate k d [20] . 
This results in the evolution equations 



d_ _d_ 
dt d dl 
d_ _ d_ 
dt a dl 



Vg(l)-V p 



v 9 (l) 



N d = -k a N d + k d N a , (la) 
N a = k a N d -k d N a . (lb) 



The right hand side of Eq. [T] describes attachment and 
detachment process. The second term on the left hand 
side accounts for the gain and loss of attached and de- 
tached polymers due to the dynamics of the polymer 
mesh, growing with velocity v g , and the polymerization 
kinetics of the filaments in the brush. The correction 
factor l/x in front of v g is due to the fact that for bent 
polymers the rigid network swallows by this amount more 
in contour length than for straight filaments. This factor 
is equal to 1 for I < x. 

Processes contributing to the growth of the rigid poly- 
mer mesh are entanglement and crosslinking of filaments 
in the brush. Both imply a vanishing v g for I — > 0, since 
short polymers do not entangle and crosslinking proteins 
are unlikley to bind to them. At the same time v g can 
not grow without bound but must saturate at some value 
^max ^ ue £ Q ra ^ e limitations for crosslinker binding. This 
suggests to take the following sigmoidal form 



v g (l)=vf ax tznh(l/T), 



(2) 



with a characteristic length scale I. 

The polymerization rate is proportional to the prob- 
ability of a gap of sufficient size d (~ 2.7 nm) between 
the polymer tip and the obstacle for insertion of an actin 
monomer |7j. This implies an exponential dependence of 
v p on the force by which the polymer pushes against 
the obstacle, 



v p (l, x) = v™* exp [-d • F d (Z, x)/k B T] . 



(3) 



Here, v™ ax ~ 500 nm s _1 [7 is the free polymerization 
velocity. For the entropic force F& we use the results ob- 
tained in Ref. [8] for D = 2, 3 spatial dimensions, where 
we take the accepted value of t v « 15 fim [21] for the 
persistence length of F-actin. 

The dynamics of the distance x between grafted end 
of the filament and the obstacle interface (see Fig. [T]) is 
given by the difference of the average v g and the velocity 
of the obstacle 

d t x = -^ J™ dlv g {l) [N a {l,t) + N d (l,t)] (4) 
dl [N a (l,t) F a (l,x)+N d (l,t) F A (l,x)] , 



(Jo 

where £ is an effective friction coefficient of the obsta- 
cle. The force F a (l,x) acting on the obstacle interface 



results from the compliance of the filaments attached 
to it by some linker protein complex, which we model 
as springs with spring constant ki and zero equilibrium 
length. This complex has a nonlinear force-extension re- 
lation which we approximate by a piece- wise linear func- 
tion; for details see the supplementary material. Let 
R\\ « l[l — 1{D — l)/4£ p ] be the equilibrium length of the 
filament. Then, the elastic response of filaments expe- 
riencing small compressional forces (x < Ru) is approxi- 
mated by a spring constant k\\ = 12k B T £ 2 /(D — [22] . 
For small pulling forces (x > i?y), the linker- filament 
complex acts like a spring with an effective constant 
&eff = hk\\/(ki + fey). In the strong force regime, the 
force-extension relation of the filament is highly nonlinear 
and diverges close to full stretching [23]. Therefore, only 
the linker will stretch out. The complete force-extension 
relation is captured by 



F a 



-k\\(x-R\\), 
-k e fi(x -R\\), 
-ki(x - I) - k e ff(l 



R\\), 



x < R\\ , 
R\\ < x < I 
x > I . 



(5) 



Finally, we specify the force-dependence of the detach- 
ment rate by 



k° d exp[-d-F a (l,x)/k B T] . 



(6) 



Here, k° d « 0.5 s _1 [12 is the detachment rate in the 
absence of forces and we have followed Ref. [20] . 
has a singularity at v p (l s ) ■ 



Eq. 



la 



x since 

the coefficient of the derivative of N d with respect to 
I is zero at l s . To illustrate the key physical features 
at that singularity, we start with the simple equation 
dtN d — di[v g (l)l/x — Vp(l,x)]N d = with x kept constant. 
Then those parts of the distribution of N d with I < l s will 
grow and catch up with l s since v g (l)l/x — v p (l, x) is posi- 
tive there, while the parts with I > l s will shorten towards 
l s . As a consequence the whole distribution will become 
concentrated at l s . To quantify this heuristic argument 
we expand v g (l)l/x — v p (l,x) up to linear order around 
l s like Vi(l — l s ) and use the method of characteristics to 
solve the equation. Starting initially with a Gaussian dis- 
tribution we obtain N d (l, t) = c(t) exp[— (I — l(t)) 2 /w(t) 2 ] 
with c(t) = coexp(^it), T(t) = l s + (lo — l 8 ) exp(— v\t) 
and w(t) = woex.p(—vit). This shows that N d evolves 
to a monodisperse distribution which is localized around 
l s . Its width decreases exponentially with time while its 
height grows exponentially. The time scale for this con- 
traction is given by [di(v g l/x — Vp)] -1 . 

Since the same kind of singularity also occurs in the 
full set of dynamic equations, Eqs. [l] we may readily 
infer that N a and N d evolve into delta-functions with 
that dynamics. This is well supported by simulations, 
and allows us to continue with the ansatz 



N d (l,t)=n d (t) S(l-l d (t)), 
N a (l,t) = n a (t) 5(l-l a (t)). 



(7a) 
(7b) 
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It defines the dynamic variables n d (t), l d (t), n a (t) and 
l a (t). Upon inserting Eqs. [7] into Eqs. [I] and Eq. [4| we 
obtain the following set of ordinary differential equations 



d t ld(t) = v p (l d ,x) 



-Vg(ld) 



d t l a (t) 

d t n a (t) 
d t x(t) 



X 

- k a — (l d 

n„, 



7 U °>n 
kd — (la 

n d 



-k d (l a ,x) n a (t) 



la), 

k a n d (t) , 



Id) , (8a) 

(8b) 
(8c) 



1 

c 

1 

N 



n a (t) F a (l a ,x)+n d (t) F d (l d ,x)] 
[v g (la) n a (t) +v g (l d ) n d (t)} , 



(8d) 



where n d (t) = N — n a (t) since we keep the total number 
of filaments fixed. 

The values of many parameters in the dynamics can be 
estimated using known properties of actin filaments. We 
choose the linker spring constant ki « 1 pN nm" 1 [12] 
and assume N = 200 [12 filaments to be crowded behind 
the obstacle. A realistic value of the drag coefficient ( is 
10 -3 pN s nm -1 but results did not change qualitatively 
for a range from 10 -5 pN s nm -1 to 1 pN s nm -1 . 

We have numerically solved Eqs. [8] in both D = 2 
and D = 3 dimensions, and found the dynamic regimes 
shown in Fig. [2j stationary states and oscillations. The 
existence of an oscillatory regime is very robust against 
changes of parameters within reasonable limits including 
the spatial dimension. We checked robustness against 
changes in the parameter values for the number of poly- 
mers N, I (see Eq.pl , ki, v™ ax and k®, in addition to the 
examples shown in Fig. [2] In general, we find that oscil- 
lations occur for v™ ax < 500 nm s _1 and within a range 
of values for k a . Note that the oscillatory region in pa- 
rameter space depends on the orientation $ of filaments 
with respect to the obstacle surface, i.e. oscillating and 
non-oscillating sub-populations of filaments may coexist 
in the same network. 

Oscillations appear with finite amplitude and period at 
the lower boundary of the oscillatory region; compare the 
example shown in Fig. [3^i. The stationary state changes 
stability slightly inside the oscillatory regime and oscil- 
lations set in with a finite period. That is compatible 
with oscillations appearing by a saddle node bifurcation 
of limit cycles. The upper boundary of the oscillatory 
region is determined by a Hopf bifurcation. An exam- 
ple of an oscillation close to that bifurcation is shown 
in Fig. [3}d. More details on the phase diagram will be 
published elsewhere [15] . 

We start with the description of oscillations in the 
phase with v g > v p: i.e., decreasing lengths x, l a and l d ; 
see Fig. [3j Then the magnitude of pulling and pushing 
forces increases due to their length-dependence. When 
the pushing force becomes too strong, an avalanche-like 
detachment of attached filaments is triggered and the ob- 
stacle jerks forward; compare the steep rise in l dl l a and x 
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FIG. 2: (color online) Phase diagram of Eqs. 8a -]8d] outlining 
stationary and oscillatory regimes with £ = 10 - pN s nm -1 
for (a-c) and (a) D = 2, & = , (b) D = 2, & = tt/4, (c) 
D = 3, & = and (d) D = 3, & = 0, C = 10" 5 pN s nm" 1 . 
/=100 nm, all other parameter values are specified in the text. 
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FIG. 3: (color online) x , l a , Id (in nm) and n a as a function 
o f tim e, as obtained from a numerical solutions of Eqs. [8a] 
with v™ ax = 300_nm s -1 and (a) k a = 0.143 s _1 (b) 
3.49 s" 1 . D = 3, Z=100 nm in both panels. 
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shown in Fig. [3] That causes a just as sudden drop of the 
pushing force. With low pushing force now, polymeriza- 
tion accelerates and increases the length of detached fila- 
ments. The restoring force of attached filaments is weak 
in this phase due to their small number. Hence, despite 
of not so strong pushing forces, the obstacle moves for- 
ward. In the meantime, some detached filaments attach 
to the surface such that the average length and number of 
attached filaments increases as well. When the detached 
filaments are long enough to notice the presence of the 
obstacle interface, they start to buckle. This, in turn, in- 
creases the pushing force and slows down the polymeriza- 
tion velocity. Therefore, the graft velocity now exceeds 
the polymerization velocity and the average lengths of 
attached and detached filaments start to decrease again 
and the cycle starts anew. The period of oscillations is 
dependent on the parameter values. It reduces from 240 s 
in Fig. [3^i to 13 s in Fig. (3)3 as k a increases from 0.143 s _1 
to 3.49 s -1 at ^ ax = 300 nm s _1 . 

The oscillations in x correspond to the saltatory mo- 
tion of the obstacle in the lab frame and the oscillations 
of its velocity since v g stays essentially constant. An il- 
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lustration is shown in Fig. [4] for a given set of parameters 
which leads to oscillations with periods of the order of 
100 s and velocity of the order of 0.7 fim s _1 . This is 
in good agreement with the results of experiments on os- 
cillatory Listeria propulsion [16]. The period of velocity 
oscillations with beads propelled by actin polymerization 
differs from those of Listeria by one order of magnitude 
(8 — 15 min [H]). Periods of that length can be obtained 
within our model upon using values for k a close to the 
lower boundary of the oscillatory regime. 
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FIG. 4: (color online) Velocity and displacement of the obsta- 
cle as a function of time with (a) k a = 0.9 s _1 , (b) k a = 1 s _1 . 
v™ ax = 750 nm s~\ < ax = 75 nm s~\ k° d = 0.1 s _ \ 
/=100 nm £ — 10 -3 pN s nm -1 and D = 3 in both panels. 



We have also studied the system when the network 
is oriented at an angle $ = 7r/4. In this case, the spring 
constant of the attached filaments parallel to n for D = 2 



reads fcjT 1 ^) = + e~ e / 2 - 1 + cos2tf(| 



4 1 12 c 



COS' 



; tf(e- 



e/2 



1) 2 ]/^B^, where e = l/l p and 



1(1 — //4^ p )cos^ [22]. For the pushing force 



|e-/») 

R\\W -- 

of a filament grafted at d = 7r/4, we use the results of 
the factorization approximation given in Ref. [8 , which is 
well valid for a stiff filament like actin. A numerical solu- 
tion of Eqs. 8a]|8d| results in the phase diagram shown in 
Fig.[2|b) with the adapted forms of and F a . The main 
effect is that one needs higher values for the attachment 
rates and lower values for v g to obtain oscillations. 

In summary, we have presented a simple and generic 
theoretical description of oscillations arising from the 
interplay of polymerization driven pushing forces and 
pulling forces due to binding of actin filaments to the 
obstacle. The physical mechanism for such oscillations 
relies on the load-dependence of the detachment rate and 
the polymerization velocity, mechanical restoring forces 
and eventually also on the cross-linkage and/or entangle- 
ment of the filament network. The oscillations are very 
robust with respect to changes in various parameters, i.e. 
are generic in this model. Therefore, complex biochem- 
ical regulatory systems supplementing the core process 
described here may rather stabilize motion and suppress 
oscillations than generate them. 

Oscillations of the velocity were observed during 
propulsion of pathogens by actin polymerization. There, 
the core mechanism described here is embedded into a 
more complex control of polymerization, which e.g. also 
comprises nucleation of new filaments and capping of ex- 



isting ones. Hence, the study presented here can not be 
expected to fully capture all features of such processes. 
Our results still agree well with respect to velocity spike 
amplitudes and periods in Listeria experiments reported 
in Refs. [16] [17]. The velocity in between spikes appears 
to be smaller in experiments than in our simulations. 
This may be accounted for in our model by including 
capping of filaments upon dissociation from the obsta- 
cle. Periods may also become longer when capping and 
nucleation were included since it would take longer to 
restore the pushing force after the avalanche like rup- 
ture of attached filaments. Altogether, qualitative and 
quantitative comparison with experiments suggests that 
our model may be a promising candidate for a robust 
mechanism of velocity oscillations in actin-based bacte- 
ria propulsion. 
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